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A  control  strategy  is  proposed  for  a  minimally-actuated  flapping-wing 
micro  air- vehicle  (FWMAV).  The  proposed  vehicle  is  similar  to  the  Har¬ 
vard  RoboFly  that  accomplished  the  first  takeoff  of  an  insect  scale  flapping 
wing  aircraft,  except  that  it  is  equipped  with  independently  actuated  wings 
and  the  vehicle  center-of-gravity  can  be  manipulated  for  control  purposes. 
Using  the  results  from  the  derivation  of  the  aerodynamic  forces  and  mo¬ 
ments  from  Part  I,  a  control  allocation  strategy  and  a  feedback  control 
law  are  designed  that  enables  the  vehicle  to  achieve  untethered,  stabilized 
flight  about  a  hover  condition.  The  control  laws  are  designed  to  make  use  of 
three  actuators,  two  of  which  control  the  angular  position  of  the  wing  in  the 
stroke  plane,  and  one  that  moves  a  bob-weight  that  manipulates  the  vehicle 
center-of-gravity.  The  Split-Cycle  Constant-Period  Frequency  Modulation 
technique,  introduced  in  Part  I,  is  used  to  allow  each  wing  to  generate  non¬ 
zero  cycle-averaged  rolling  and  yawing  moments.  The  technique  achieves 
this  objective  by  varying  the  frequency  of  the  oscillators,  that  drive  each 
wing  throughout  the  wing-beat  cycles,  such  that  the  dynamic  pressure  act¬ 
ing  on  each  wing  during  the  upstroke  is  different  from  that  which  acts  on  the 
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wing  during  the  downstroke.  Pitching  moment  and  high  speed  translation 
are  controlled  by  varying  the  vehicle  center-of-gravity  using  a  bob-weight 
actuator. 


I.  Introduction 

In  Part  I1  of  this  series,  expressions  for  the  instantaneous  and  cycle  averaged  aerody¬ 
namic  forces  and  moments  of  a  flapping  wing  biomimetic  vehicle,  with  features  similar  to 
the  Harvard  Robofly  described  by  Wood,2  were  derived.  The  proposed  vehicle  makes  use 
of  three  actuators,  two  of  which  are  used  to  independently  vary  the  position  of  the  wing 
spars  in  the  stroke  plane,  while  the  third  is  used  to  control  the  position  of  a  bob-weight. 
The  planforms  are  free  to  passively  rotate  about  the  spar  within  limits  prescribed  by  a  spe¬ 
cially  designed  hinge,  the  details  of  which,  are  provided  in  Part  I  of  this  series.  Because 
of  the  limited  number  of  actuators  and  the  periodic  nature  of  the  aerodynamic  forces  and 
moments,  a  control  strategy  has  been  devised  that  will  allow  5  degree  of  freedom  control  of 
the  fuselage  to  be  achieved.  The  split-cycle  constant-period  frequency  modulation  method 
that  enables  this  motion  is  a  method  which  allows  the  frequency  of  the  oscillator  that  con¬ 
trols  the  instantaneous  wing  position  to  be  varied  at  the  endpoints  of  each  stroke,  in  such 
a  way  as  to  preserve  the  period  of  a  temporally  symmetric  wing  beat.  Thus,  the  angular 
velocity  and  dynamic  pressure  can  be  varied  throughout  each  cycle.  This  feature  enables 
asymmetric  cycle  averaged  forces  and  moments  to  be  generated  that  can  be  used  for  maneu¬ 
vering  and  control.  There  are  two  input  parameters  to  each  wing  beat  waveform  generator. 
These  parameters  are  called  the  fundamental  wing  beat  frequencies,  ulw  and  corw,  and  the 
split-cycle  parameters,  Slw  and  5rw-  The  fundamental  frequency  controls  the  frequency  of 
the  temporally  symmetric  portions  of  the  wing  beat,  while  the  split-cycle  parameters  control 
the  temporally  asymmetric  portions  of  the  wing  beat. 

In  the  present  paper,  the  analysis  that  was  started  in  Part  I  is  continued  by  describing 
the  split-cycle  formulation  and  analyzing  the  effects  of  a  bob-weight  upon  the  body-axis 
forces  and  moments.  Expressions  for  the  control  derivatives  of  the  cycle-averaged  forces  and 
moments  associated  with  variations  in  the  fundamental  wing  beat  frequencies,  split-cycle 
parameters,  and  bob-weight  position  are  derived.  The  expressions  are  linearized  about  a 
hover  flight  condition  and  a  control  allocation  strategy  is  proposed.  A  flight  control  law  is 
formulated  based  on  the  linearized  cycle-averaged  model,  and  is  implemented  on  a  dynamic 
model  that  includes  instantaneous  variations  in  the  forces  and  moments  to  verify  the  efficacy 
of  the  cycle-average  based  control  design.  Areas  where  dynamics  and  control  can  influence 
the  design  of  the  aircraft  are  highlighted  throughout  the  paper. 
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II.  Equations  of  Motion 


A  standard  set  of  6  degree-of- freedom  rigid  body  equations  of  motion  are  employed.3 
These  are  the  same  equations  that  describe  the  motion  of  a  rigid  body  aircraft;  however, 
because  the  angle  of  each  wing  in  the  stroke  plane  is  a  periodic  function  of  time,  the  forces 
and  moments  that  drive  the  equations  consist  of  multiple,  independent,  and  variable  periodic 
functions  of  time.  We  assume  that  the  wings  are  massless;  therefore,  there  exists  no  inertial 
coupling  between  the  fuselage  and  the  wings.  We  expect  that  the  inertial  coupling  between 
the  wings  and  body  will  be  small  because  in  the  Harvard  Robofly,  the  wings  constitute  less 
than  1%  of  the  total  vehicle  weight.2  The  equations  of  motion  are 
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where  I  is  the  inertia  matrix,  u>  =  \p  q  r]  is  the  angular  rate  vector,  L ,  M,  N  are  the  rolling, 
pitching,  and  yawing  moments,  u,  v,  w  are  the  translation  velocities,  X,  Y,  Z  are  the  body 
axes  forces,  m  is  the  vehicle  mass,  g  is  the  acceleration  of  gravity,  x,  y,  z  are  the  vehicle 
positions  with  respect  to  an  inertial  frame,  and  q  =  9o  +  9Y  +  9‘ij  +  93 k  is  a  quaternion 
used  to  perform  the  3-2-1  standard  Euler  transformation.  Also,  =  HfT  are  rotation 
matrices  that  transform  from  body  to  inertial  axes  and  back.  The  equations  of  motion  are 
driven  by  the  instantaneous  forces  and  moments  presented  in  Part  I1  of  this  series.  These 
blade  element  based  aerodynamic  forces  and  moments  do  not  account  for  leading-edge  vortex 
effects  and  the  effects  of  wake-capture  that  have  been  observed  in  experiments  and  therefore, 
in  practice,  there  exists  a  significant  amount  of  uncertainty  that  must  be  accommodated  by 
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any  feedback  controller  designed  based  on  blade-element  results. 


III.  Split-Cycle  Fundamentals 

Varying  the  fundamental  wing  beat  frequency  allows  control  over  the  translational  and 
rotational  degrees  of  freedom  of  the  vehicle.  When  the  split-cycle  parameters  5lw  and  8rw 
are  zero,  the  left  and  right  wings  beat  symmetrically  and  multi-axis  control  of  the  vehicle  is 
not  possible.  This  section  describes  the  split-cycle  parameters  and  shows  how  the  wing  beat 
frequency  is  modified  for  the  upstroke  and  downstroke. 

The  fundamental  idea  is  to  piece  together  two  cosine  waves  of  different  frequencies  over 
one  complete  wing  beat  cycle.  Either  the  upstroke  or  the  downstroke  will  be  impeded,  while 
the  opposite  stoke  is  advanced.  The  development  of  the  split-cycle  parameters  for  the  right 
wing  is  exactly  the  same  as  for  the  left  wing.  Therefore,  the  subscripts  for  right  wing  and 
left  wing  will  be  eliminated  in  the  following  derivations. 

For  the  first  half  of  the  cycle,  the  frequency  is  altered  from  cu  to  u  —  8.  Thus,  the  position 
of  the  wing  during  the  upstroke  is  defined  by 

4>u  =  cos  [(a;  —  5)  t\  (5) 

where  5  is  the  split-cycle  parameter  and  c o  is  the  symmetric  or  fundamental  wing  beat 
frequency.  For  8  >  0,  this  would  produce  an  impeded  or  lower  frequency  upstroke,  as 
compared  to  the  fundamental  wing  beat  frequency.  For  the  downstroke,  the  position  of  the 
wing  becomes 

4>d  =  cos  [(cu  +  a)  t]  (6) 

where  a  is  the  frequency  modification  parameter,  which  is  dependent  upon  5.  In  order  to 
compute  the  fundamental  frequency  change  for  the  second  half  of  the  cycle,  recall  that  the 
combination  of  the  two  altered  frequency  strokes  must  have  the  same  period  as  the  unaltered 
frequency  stroke.  Then,  for  an  upstroke  frequency  of  uj  —  8  and  a  downstroke  frequency  of 
uj  +  cr,  the  requirement  becomes 

Tv  +  Td  =  T  (7) 

where  Tjj  =  -jT  is  the  period  of  the  upstroke,  To  —  is  the  period  of  the  downstroke,  and 
T  =  —  is  the  full  cycle  period.  Since  the  upstroke  and  downstroke  only  occur  for  one-half 
of  their  respective  periods,  Equation  7  becomes 

1  /  2vr  \  1  /  2vr  \  2tt 

2  (yvy)  +  2  (yrry)  ~~  (8) 
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Solving  for  o  yields 


Su> 

uj  —  25 


(9) 


Equation  9  shows  that  a  is  a  function  of  the  symmetric  frequency,  uj,  and  the  split-cycle 
parameter,  5.  Thus,  split-cycle  frequency  modulation  is  fully  defined  by  the  parameter  5. 
There  are  two  equations  that  provide  limits  on  <5,  namely, 


uj  —  <5  >  0  cu  T  ex  >  0 


(10) 


Equation  10  ensures  that  the  wing  beat  frequencies  are  always  positive.  The  first  part  reduces 
to  S  <  tv.  For  the  second  part,  substitute  Equation  9  into  Equation  10  and  simplify  to  obtain 
the  following  condition 


uj  (uj  —  5) 


>  0 


in) 


u:  —  25 

Three  cases  occur  here.  Note  that  the  numerator  of  Equation  11  is  always  positive  since 
5  <  uj  and  uj  >  0.  Consider  the  denominator  of  Equation  11.  Case  1  is  for  5  >  j. 


uj  —  25  <  0  — >  uj  +  a  <  0 


(12) 


Obviously,  the  frequency  of  the  downstroke  cannot  be  less  than  zero  and,  therefore,  5  cannot 
be  larger  than  j.  Case  2  is  for  5  —  j.  Then, 

5  —  —  — >  uj  —  25  =  0  — ■>  uj  +  a  —  oo  (13) 


As  with  case  1,  case  2  yields  an  unrealistic  result.  Case  3  is  for  <5  <  j.  Then, 

<5  <  —  — >  uj  —  25  >  0  — »■  uj  +  a  >  0 


(14) 


For  case  3,  the  frequency  of  the  downstroke  is  larger  than  zero,  as  required.  Since  this 
requirement,  namely,  5  <  is  more  restrictive  than  that  obtained  from  the  first  part  of 
Equation  10,  we  conclude  that  5  must  satisy 


<5  <  |  (15) 

If  5  =  0,  the  fundamental  wing  beat  frequency  for  both  the  upstroke  and  downstroke  is 
achieved. 

It  should  be  noted  that  5  <  0  is  also  an  option.  In  this  case,  the  upstroke  is  advanced 
and  the  downstroke  is  impeded  (just  the  opposite  of  the  <5  >  0  case).  There  is  no  lower  limit 
on  5,  which  can  again  be  determined  using  Equation  10.  Clearly,  for  5  <  0,  the  first  part  of 
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Equation  10  is  satisfied.  Consider  the  following  for  the  oo  +  &  >  0  requirement: 


Then, 


800  800  —oo 

Inn  a  =  lim  - —  =  Inn  — —  =  — — 

<5— >— oo  <5— >— oo  00  —  2 8  <5— >— oo  — 2o  2 


00  00 

lim  oo  +  a  =  oo  —  —  =  — 
5— >— oo  2  2 


(16) 


(17) 


and  thus,  the  downstroke  frequency  is  greater  than  zero  for  any  8  <  0.  Therefore,  the  only 
active  constraint  for  the  split-cycle  is  <5  <  j. 

One  other  requirement  needs  to  be  taken  into  account.  Equation  6  must  to  be  modified 
to  ensure  that  the  split-cycle  wing  position  is  continuous.  In  other  words,  the  wing  position 
at  the  end  of  the  upstroke  must  equal  the  wing  position  at  the  beginning  of  the  downstroke. 
The  upstroke  ends  at  time  t  =  .  Substituting  this  time  into  Equation  5  yields 


0{/|t=  -  =  COS 


(w-5) 


7T 


oo  —  8 


=  COS  7T  =  — 1 


(18) 


Evaluating  Equation  6  at  this  time  produces 


<f)D\t=*  =  cos 


(oo  +  a) 


71 


oo  —  8 


=  cos 
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IT 
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c* 
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=  COS 
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On 

oo  —  28 

(19) 


Considering  the  argument  of  the  cosine  function  and  setting  it  equal  to  7 r  (the  result  of 
Equation  18)  gives 

U7rnr  =  n  -»•  28n  =  0  -»•  5  =  0  (20) 

lo  —  28 

Hence,  there  will  be  a  mismatch  of  wing  positions  at  the  transition  from  the  upstroke  to  the 
downstroke  if  8  ^  0.  To  eliminate  this,  the  downstroke  can  be  phase  shifted.  The  form  of 
the  downstroke  wing  position  function  is  altered  from  Equation  6  to 


4>d  =  cos  [(a;  +  a)  t  +  £] 

Evaluating  Equation  21  at  time  t  =  produces 

oo  (u  —  8)  f  71 


(21) 


(j)D\t=  =  COS 


(oo  +  a) 


7 r 

lo  —  8 


+  £ 


=  cos 


oo  —  28  V  c o  —  8 


0071 

=  cos 

- ^7 

oo  —  28 

(2 

Setting  the  argument  of  the  cosine  in  Equation  22  equal  to  7r  produces 

-27r5 


0071 


oo  —  28 


f  =  TT  ->  i  = 


oo  —  28 


(23) 


Utilizing  the  phase  shift  in  Equation  23  guarantees  that  the  wing  position  at  the  end  of  the 
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upstroke  is  equal  to  the  wing  position  at  the  beginning  of  the  downstroke. 

In  summary,  the  upstroke  and  downstroke  wing  positions  are  given  by 

( i>u  =  cos  [(cu  —  5)  t]  (24) 

c t>D  =  cos  [(a;  +  a)  t  +  £]  (25) 

where  a  =  and  £  =  -iff-  The  requirement  on  the  split-cycle  parameter  is 

<5  <  f  (26) 

and  6  can  take  on  any  negative  value. 

Figures  1  and  2  show  wing  beat  patterns  for  S  >  0  and  6  <  0  along  with  a  symmetric 
wing  beat.  In  Figure  1,  the  upstroke  is  impeded,  while  the  downstroke  is  advanced,  while  in 
Figure  2,  the  upstroke  is  advanced  and  the  downstroke  is  impeded.  The  split-cycle  parameter, 
<5,  along  with  u  for  the  right  wing  and  left  wing  become  control  parameters  in  the  control  law. 
Hence,  the  wing  position  control  parameters  are  u>rw,  w lw >  $rw,  and  &lw-  For  utilization 
in  a  control  law,  the  control  derivatives  associated  with  the  fundamental  and  split-cycle 
frequencies  must  be  evaluated.  The  next  section  addresses  this  task. 

1 

0.8 
0.6 
0.4 
0.2 
I  0 

-e- 

-0.2 
-0.4 
-0.6 
-0.8 
-1 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Normalized  Time 


5  >  0 


Figure  1.  Split-cycle  results  for  S  >  0. 
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5  <  0 


Figure  2.  Split-cycle  results  for  S  <  0. 


IV.  Control  Derivatives 

The  parameters  used  to  control  the  aerodynamic  forces  and  moments  are  the  fundamental 
wing  beat  frequencies,  corw  and  ujlw,  and  the  split-cycle  parameters,  5rw  and  Slw-  For  con¬ 
trollability  analysis  and  control  synthesis,  the  sensitivity  of  each  cycle-averaged  generalized 
force  to  each  control  input  parameter  must  be  determined. 


A.  General  Aerodynamic  Control  Derivative  Expressions 

The  partial  derivatives  of  each  cycle-averaged  force  and  moment,  derived  in  Part-1,  with 
respect  to  each  wing  beat  control  input  parameter  are  provided  below 

_ q 

dFxRW  _  rw5 rw{u rw  —  Srw) 

OSrw  (urw  ~  2 Srw)2 


dFXRw  _  ^l(^rw  ^rw^rw  +  ^mrw&rw  $rw) 
Sujrw  i^Rw  ~  2 5rw)2 


(28) 
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dMyBLW  _  Ji(1)ulw  [^LX^pp  cos  a  +  kDx™p  sin  a  +  kDAx f  ]  {282LW  -  2 8lwojlw  +  u 


2  1 
LW) 


9 SLw  (^lw  ~  2 8lw)2 

ki,Azp  8 i,w^ lw(,^ lw  ~  8lw ) 

—  28LW)2 


(42) 


dMyLW  _Ji(1)8lw  [kLx™p  cosa  + kDxcp 


WP 


du 


lvlw 

LW 


sin  a  +  kDAxp\  (: 282LW  -  4 8lwdLw  +  u 2LW ) 


{x>lw  —  2 5lw ) 


kxAzp  (8\w  4  8‘fjyyU>Lw  T  ^lw^Ilw  ^ 


lw  ) 


{x>lw  —  28lw) 


(43) 


d Mz^w  _  — 2A;L(y^pJi(l)  +  ^)drw8rw(urw  —  <5RM/) 

<9<5rw  (^rw  —  28  rw)2 

9MzRW  _  —  2ki(y^p Ji(l)  +  ^)(^Rw  —  4 8RwuR\y  +  4u;rwA^w  —  ^w) 
duRw  (uRw  ~  28 Rw)2 


dMzfw  _  2kL{y^pp JiiX)  +  f  )ulw8lw(ulw  —  8lw) 

95lw  {x>lw  —  2  8lw)2 

9MZ mv  _  2^i(^cpP^i(l)  +  f  )(ujlw  ~  4 8lw^lw  +  4a )lw8\w  ~  $lw) 

dioLW  {uLW  —  2  8lw)2 


(44) 

(45) 

(46) 

(47) 


B.  Aerodynamic  Control  Derivatives  about  Hover 

The  control  of  this  vehicle  in  the  vicinity  of  hover  is  of  considerable  interest.  Therefore,  the 
control  derivatives  at  the  hover  condition  where  c oRw  =  dlw  =  o;G  and  8Rw  =  8lw  =  0 
are  evaluated,  where  ua  =  {)  •  Of  course,  this  condition  assumes  that  the  nominal 

center-of-gravity  of  the  vehicle  and  wing  root  hinges  are  aligned  such  that  A zR  =  Azf  =  0. 
The  control  derivatives  about  hover  are  presented  in  Table  1. 
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Control  Derivative  ^\u=w0,s= o 

9FXRW 

96rw 

=  0 

urw=u0,Srw= 0 

9Fxrw 

9urw 

=  kLua 

urw=u0,Srw= 0 

9FXLw 
96  lw 

=  0 

L  W  — 1 Mo  , s  L  W  =' 0 

9Fxrw 

9ujlw 

=  kLu0 

ut.w—u0,6t.w—9 

rw  _  RW  _  dFVLW  _  ®FyLW  _  g 

dujRvv  dS  r\a/  dto  r,w  OSt.w 

9FZ  rw 
96rw 

—  —  kjj  Ji(l)u0 

ojrw=Vo,6rw=0 

9FZ  rw 
9urw 

=  0 

u>rw=Wo,6rw=0 

9Fzlw 
96  lw 

=  —kr>Ji(l)uj0 

ut,w=u0,6t,w= 0 

9Fzlw 

9u>lw 

=  0 

ur,w=u0,6r,w= 0 

9Mxrw 

96rw 

=  -\kDu0(yWp  +  wJi(l)) 

urw—  Vo,6rw—' 0 

9Mxrw 

9<jjrw 

=  0 

urw=Vo,6rw=0 

dMx  lw 
96lw 

=  \kDu0(yWp  +  wJi(l)) 

vlw—v0,6lw— 0 

dMx  lw 
9ujlw 

=  0 

u>r,w=‘x’o,6r,w=0 

9M~yRW 

96rw 

=  Ji(l)u!0  [kLX™pP  cos  a  +  /cD(x^p  sin  a  +  Axj|)] 

urw=u0,Srw= 0 

9M-vBrw 

9ujrw 

=  kL£±ZRU0  =  0  for  A zp  =  0 

IjJRW=(jJo,6pw=  0 

9M-ylw 
96  lw 

=  Ji(l)cuD  \kLx^p  cosa  +  kD(x^pp  sina  +  Axf)] 

u)r,w=‘x>o,6r,w=0 

9M~ylw 

9ll)Lw 

=  kLAzpuQ  =  0  for  A zp  =  0 

<jJlw=<jJo,6t,w= 0 

9MzrW 

96rw 

=  0 

u>rw=u>o,6rw=0 

9MzrW 

9lorW 

=  -2 kLu,„  [y«/M  1)  +  f] 

WBW=Wo,ORW=0 

9MZLw 
96  lw 

=  0 

u>t,w=u>o,6t,w=Q 

9MZLw 

9ulw 

=  2*Wc,[<pJ1(l)  +  f] 

uu,w=w„,ar,w=0 

Table  1.  Aerodynamic  control  derivatives  about  hover. 
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C.  Bob- weight  Control  Derivatives 


It  is  important  to  note  that  variations  in  symmetric  and  split  cycle  frequency  parameters 
are  not  effective  in  generating  independent  changes  to  the  pitching  moment.  Using  only 
the  aerodynamic  control  input  parameters,  changes  in  pitching  moment  are  accompanied  by 
undesirable  changes  in  other  forces  and  moments.  It  is  for  this  reason  that  the  bob-weight 
actuator  is  introduced.  The  bob-weight  is  taken  to  lie  in  the  vehicle  plane  of  symmetry  and 
at  zero  deflection  is  coincident  with  the  center-of-gravity  of  the  vehicle  without  the  bob- 
weight.  The  bob- weight  need  not  be  dead- weight.  Any  necessary  vehicle  component  that  is 
not  sensitive  to  movement  could  serve  the  purpose,  e.g.,  battery,  power  converter,  etc.  The 
weight  would  be  mounted  near  the  tip  of  a  bimorph  piezoelectric  actuator.  The  tip  deflection 
is  designed  to  primarily  consist  of  translation  in  the  z-body  axis  direction,  while  translation 
in  the  x-body  direction  would  be  so  small  as  to  be  negligible.  Hence,  the  bob-weight  motion 
is  constrained  to  be  an  arc  in  the  Xb,  Zb  plane. 

Let  m0  be  the  total  vehicle  mass  exclusive  of  the  bob-weight  and  bob-weight  actuator 
and  xfgo,0,z^go  be  the  location  of  the  center  of  mass  in  body  coordinates.  Let  mBw  denote 
the  mass  of  the  bob-weight  and  actuator  assembly  and  x^w,  0,  z^w  be  the  location  of  the 
center  of  mass  of  this  assembly  in  body  coordinates.  Because  motion  is  constrained  to  an 
arc  in  the  Xb,  Zb  plane,  the  y-body  locations  of  these  centers-of-mass  are  zero,  that  is,  no  out 
of  plane  motion  occurs  and  the  radius  of  the  arc  at  the  tip  of  the  actuator  is  large  relative  to 
the  tip  displacement.  Because  of  the  nature  of  the  tip  motion  of  the  bimorph  actuator,  it  is 
assumed  that  motion  of  the  bob-weight  does  not  alter  the  x-location  of  the  center  of  mass  of 
the  bob-weight  assembly,  hence,  XpjW  is  constant.  Defining  x^g,  0,  as  the  location  of  the 
center  of  mass  of  the  entire  vehicle,  including  the  bob-weight,  this  assumption  implies  that 
x^y  is  fixed  for  a  given  bob-weight  mass.  Wood2  indicates  that  the  mass  of  each  wing  is  less 
than  1%  of  the  gross  vehicle  weight;  therefore,  it  is  also  assumed  that  the  mass  of  each  wing 
is  negligible.  Under  the  above  conditions,  the  x-  and  z-body  locations  of  the  center  of  mass 
are  given  by 


Xcg0mo  +  x BWmBW 
m0  +  rriBw 

Zcg0mo  +  zBWmBW 

m0  +  mBW 


(48) 

(49) 


Again,  for  a  fixed  bob-weight  mass,  the  x  location  of  the  center  of  mass  is  independent  of 
the  movement  of  the  bob-weight.  It  is  important  to  note  that  the  aerodynamic  moments  are 
taken  about  the  body-frame  and  are  independent  of  the  center-of-gravity  location,  even  when 
the  center-of-gravity  is  not  coincident  with  the  origin  of  the  body-axes  system.  Therefore, 
the  moments  about  the  body-axes  are  comprised  of  aerodynamic  moments  and  moments 
due  to  center-of-gravity  offset.  In  other  words,  the  total  cycle-averaged  moment  about  the 
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body-axes  coordinate  system,  M  ,  become 


ME  =  ME  ro  +  M®  (50) 

_ g  _ g 

where  Maero  is  the  cycle-averaged  moment  produced  by  the  wings  and  Mcg  is  the  cycle- 
averaged  moment  due  to  center-of-gravity  offset  from  the  origin  of  the  body-axes.  The 
cycle-averaged  moment  due  to  center-of-gravity  offset  is 


x  F®  dt 


(51) 


X 


B 


0  z: 


where  the  integral  of  each  component  of  the  integrand  is  evaluated  and  r^g.  = 
is  the  position  vector  from  the  center-of-gravity  to  the  origin  of  the  body-axes  system  anc 
xB ,  zcn  are  provided  in  Equations  48  and  49.  The  force  is  gravitational  and  can  be  written 


Jcgi  “eg 

as 


t^B  _  p>B 

cg  — 


~(m0  +  mBw)g  0  0 


(52) 


where  Rf5  is  a  rotation  matrix  from  the  inertial  frame  to  the  body  frame  and  g  is  the  accel- 


Fb  Fb  F, 


-*B 


Ucg 


eration  due  to  gravity.  Let  this  force,  in  general,  be  represented  as  = 

Then,  by  performing  the  cross  product  operation,  the  moment  due  to  center-of-gravity  offset 

is 

r 

— 


— B  C 0 


'0 


~B  pB 


eg  Vcg 

pB  _  B  pB 
“cg±  xcg  u'cg±  zcg 

B  pB 

C9  Vcg 


dt 


(53) 


Assume  that  the  position  vector  and  force  in  Equation  52  are  constant  over  a  full  wing  beat 
cycle.  In  light  of  these  assumptions,  the  cycle-averaged  moments,  computed  by  evaluating 
the  integrals  in  Equation  53  are 


_„B  pB 

eg  Vcg 

J  pB  _  B  pB 

“cg±  Xcg  eg1  Zcg 

B  pB 

eg  Vcg 


(54) 


It  is  assumed  that  a  linear  relationship  exists  between  the  voltage  applied  to  the  bob-weight 
actuator,  Vbw,  and  the  bob- weight  tip  deflection.  Therefore, 


dzB\v 

dVsw 


ksw 


(55) 


Additionally,  since  it  is  assumed  that  x^w  is  a  constant,  =  0.  Then,  the  control 
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effectiveness  of  the  bob-weight  actuator  can  be  expressed  as 

<9Mcg  _  <9Mcg  Ozbw 
OVbw  dzgW  9Vrw 


(56) 


Substituting  Equation  49  into  Equation  53  and  computing  the  derivatives  as  specified  in 
Equation  56  yields 


0M° 

OVbw 


(  _z 


mBW 


rpB  I? 

in0+mBw 


\ 


Fb  k 

m0+mBW  Xcg 

0 


BW 


(57) 


Note  that  when  the  body  axis  system  is  aligned  with  the  inertial  axis  system,  the  attitudes 
of  the  vehicle  are  zero  and  Rf  =  I.  Thus,  FBcg  =  =  0  and  FBcg  =  —(m0  +  msw)g- 

As  shown  in  the  prior  section,  there  are  control  effectiveness  parameters  for  the  change  in 
moments  due  to  fundamental  wing  beat  frequency  and  split-cycle  parameters.  Equation  57 
shows  that  additionally,  there  are  control  derivatives  for  changes  in  rolling  and  pitching 
moments  due  to  changes  in  voltage  applied  to  the  bob-weight  actuator.  These  expressions 
will  be  used  in  construction  of  the  control  effectiveness  matrix  in  the  next  section. 


V.  Control  Allocation 

Form  the  control  effectiveness  matrix,  whose  elements  consist  of  the  aerodynamic  control 
derivatives  from  Table  1  and  the  bob-weight  control  derivatives  from  Equations  57,  such  that 


"a  fxb~ 

p  1 

A  TyB 

$RW 

AEf 

Urw 

a  -n~B 

A  Mx 

=  ba 

8lw 

\  r  B 

A  My 

*  -t-z-B 

ulw 

Vbw 

[amz  \ 

where 


0 

Bl2 

0 

#14 

0 

0 

0 

0 

0 

0 

kDJ\(l)u0 

0 

— ko  Jl(l)cdo 

0 

0 

(59) 

’o{y™p  +wJi(i)) 

0 

0 

B45 

Bsi 

0 

B53 

0 

B55 

0 

B&2 

0 

#64 

0 
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where  B12  =  Bu  =  kLu0,  B6 2  =  -BM  =  -2kLuj0  [y™p  J^l)  +  f], 

B51  =  Ji(lK  [&l^P  cos  a  +  kD(x^p  sin  a  + A^g)] .  Additionally, 

^53  =  ^i(l  Vo  ^r^Pcos«  +  kn{x^pp  sin  a  +Axf )] ,  B45  =  -kBW  (m'on+BJBW)  Fycg  and 

-B55  =  k bw  ( — '+ B,nBW  )  Fx,,„  •  I11  order  to  compute  the  control  vector,  a  pseudo-inverse  solution 
to  Equation  58  is  utilized. 

Examination  of  this  matrix  reveals  what  is  possible  to  achieve  by  manipulating  the  four 
aerodynamic  control  input  parameters  and  the  bob- weight.  The  rank  of  this  matrix  is  five. 
Most  obviously,  it  is  not  possible  to  generate  direct  side  forces  on  the  vehicle  using  the 
five  control  parameters.  The  next  observation  is  that  rows  associated  with  the  z-force  and 
pitching  moment  would  be  linearly  dependent  for  A zp  =  A zp  =  0  and  =  Axp  if  it 
were  not  for  the  presence  of  the  bob-weight  actuator.  These  conditions  would  exist  on  a  well 
designed  vehicle  because  of  the  plane  of  symmetry,  and  the  fact  that  from  a  vehicle  design 
perspective,  it  is  desirable  that  the  z-offset  between  the  nominal  center-of-gravity  and  the 
wing  root  hinge  points  be  zero.  This  arrangement  eliminates  a  source  of  undesirable  non- zero 
pitching  moment  associated  with  uQ.  Such  pitching  moments  would  tend  to  force  the  vehicle 
out  of  trim  and  destabilize  the  hover  condition  unless  active  control  were  applied.  The  split- 
cycle  control  concept  was  conceived  to  provide  a  set  of  control  inputs  that  affect  the  vehicle 
forces  and  moments  in  one  axis  at-a-timc  in  order  provide  a  high  level  of  decoupling.  The 
sparsity  of  the  control  effectiveness  matrix  shows  that  the  concept  goes  far  towards  achieving 
this  desired  effect.  Notably: 

•  Symmetric  manipulation  of  the  fundamental  wing  beat  frequencies  provides  x-body 
axis  control. 

•  Manipulation  of  the  split-cycle  parameters  has  no  effect  upon  the  x-body  axis.  Physi¬ 
cally,  this  is  because  the  increase  in  dynamic  pressure  on  one  half  of  the  stroke  is  offset 
by  a  decrease  in  dynamic  pressure  on  the  other  half,  resulting  in  zero  net  change  in 
body  force  when  averaged  over  a  cycle. 

•  Zero  z-body  force  is  achieved  when  the  split-cycle  parameters  are  zero,  which  is  a 
result  of  drag  cancelation  on  the  up  and  down  stroke  when  the  wing-beat  cycle  is 
temporally  symmetric.  Non-zero  z-body  forces  are  generated  when  the  strokes  are 
temporally  asymmetric  across  the  wing  beat  cycle.  When  the  split-cycle  parameters 
are  varied  symmetrically,  they  can  be  used  to  cause  the  vehicle  to  translate  forwards 
or  backwards. 

•  Rolling  moments  can  be  generated  by  changing  the  split-cycle  parameters  and  are 
independent  of  the  other  control  inputs.  The  split-cycle  parameters  create  a  non- zero 
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cycle  averaged  drag  on  each  wing  due  to  temporal  asymmetry  over  the  wing  beat  cycle. 
This  non- zero  drag  acting  over  a  moment  arm  creates  the  rolling  moment. 

•  As  noted  above,  pure  z-body  forces  can  be  generated  by  symmetric  variation  of  the 
split-cycle  parameters,  which  does  not  introduce  a  rolling  moment  due  to  the  opposite 
signs  and  equal  magnitude  of  the  rolling  moment  split-cycle  control  derivatives. 

•  Pitching  moments  are  generated  by  manipulation  of  the  split-cycle  parameters  and 
will  be  linearly  correlated  with  the  z-force  used  to  control  fore/aft  translation.  The 
pitching  moment  that  is  generated  tends  to  tilt  the  stroke  plane  in  the  direction  of 
translation;  thereby  leading  to  an  increase  in  fore  or  aft  velocity. 

•  Yawing  moments  can  be  generated  by  asymmetric  variations  in  the  fundamental  wing 
beat  frequencies.  Symmetric  variations  cause  a  change  in  x-body  force;  thus,  by  ap¬ 
propriate  control  allocation,  the  x-body  force  and  yawing  moment  can  be  decoupled. 
Furthermore,  yawing  moment  is  uncoupled  from  all  other  forces  and  moments,  when 
A Zft  =  A zf  =  0,  which  should  be  the  case  by  design. 

Finally,  note  that  the  rows  associated  with  the  rolling  and  yawing  moments  are  linearly 
independent  from  all  other  rows,  meaning  that  they  can  be  independently  controlled  using 
the  four  aerodynamic  control  inputs.  In  short,  with  the  four  aerodynamic  controls,  the  x- 
body  and  z-body  forces  and  rolling  and  yawing  moments  can  be  manipulated.  For  high  speed 
translation  longitudinally,  a  pitch-to-translate  strategy  can  utilized  using  the  bob- weight. 

VI.  Control  Law  Implementation 

Utilizing  the  analysis  presented,  a  simulation  of  the  flapping  wing  micro  air  vehicle  was 
developed.  Controllers  were  designed  to  track  altitude,  forward  translation,  roll,  pitch,  and 
yaw  commands.  The  choice  of  a  split-cycle  frequency  modulated  control  strategy  led  to  a 
nearly  decoupled  cycle-averaged  control  input  matrix.  Because  of  the  structure  of  the  control 
input  matrix,  and  the  fact  that  the  aircraft  near  hover  is  effectively  described  by  a  bank  of 
decoupled  integrators,  it  is  possible  to  use  simple  SISO  methods  to  design  feedback  control 
loops  for  each  controllable  axis.  A  block  diagram  of  the  control  structure  is  shown  in  Figure  3. 
Five  axes  are  actively  controlled  or  regulated:  altitude,  forward  position,  roll  attitude,  pitch 
attitude,  yaw  attitude,  and  their  derivatives.  Each  axis  is  tuned  to  achieve  a  second-order 
response  using  a  combination  of  rate  and  position  feedback  and  reference  commands  to 
specify  desired  accelerations.  These  acceleration  commands  are  multiplied  by  the  mass  (in 
the  case  of  generating  a  desired  force)  or  the  appropriate  moment  of  inertia  term  (in  the 
case  of  generating  a  desired  moment).  The  SISO  controllers  are  designed  as  though  the 
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desired  forces  and  moments  could  be  generated  as  continuous  and  steady  functions  of  time; 
however,  the  physics  of  the  flapping  wings  dictates  that  the  forces  will  be  periodic  functions 
of  time.  The  matrix  Ba  describes  the  relationship  between  the  cycle-averaged  forces  and 
moments  and  the  fundamental  frequency,  split  cycle  parameters  and  bob-weight  position. 
A  pseudo- inverse  control  allocation  scheme  then  determines  the  necessary  control  inputs  to 


Figure  3.  Five  degree-of-freedom  split-cycle  controller. 

produce  the  desired  cycle-averaged  forces  and  moments.  The  control  allocator  computes 
the  fundamental  frequencies  and  split-cycle  parameters  based  on  the  instantaneous  force 
and  moment  commands.  The  fundamental  frequencies  and  split-cycle  parameters  govern 
the  behavior  of  oscillators  that  generate  time  varying  wing  position  commands.  Inputs  to 
the  oscillators  are  held  constant  over  the  course  of  a  wing  beat  in  order  to  avoid  changing 
the  shape  of  the  split-cycle  waveform.  The  logic  used  to  accomplish  this  function  is  called 
a  cycle  zero-order  hold  (CZOH).  The  output  of  the  oscillators  are  applied  to  the  bimorph 
piezoelectric  actuators  and  it  is  assumed  that  the  actuator  bandwidth  is  sufficient  to  track  the 
wing  position  commands.  The  instantaneous  aerodynamic  forces  and  moments  are  computed 
using  a  modified  blade-element  model  and  applied  to  the  6  DOF  equations-of-motion.  The 
subsequent  section  shows  the  results  of  applying  the  split-cycle  control  law  to  a  simulation 
of  an  aircraft  that  is  driven  by  the  instantaneous  forces  and  moments  produced  by  the 
blade-element  aerodynamic  model. 


VII.  Results 


Utilizing  the  analysis  presented  here  and  in  Part  I,1  a  simulation  of  the  flapping  wing 
micro  air  vehicle  was  produced.  Control  laws  were  designed  for  altitude,  forward  translation, 
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rolling,  pitching,  and  yawing  moments.  Figure  4  shows  the  three-dimensional  trajectory 
that  the  vehicle  is  commanded  to  follow.  The  vehicle  begins  at  (0,0,0),  rotates  to  align 
heading  with  point  1  and  translates  to  point  1,  then  performs  similar  heading  alignment  and 
translation  to  move  to  points  2,  3,  4,  and  finally  5.  Some  of  the  vehicle  properties  are  shown 

in  Table  VII. 


Variable 

Value 

LI  nits 

rn0 

60 

mg 

Tflbw 

20 

mg 

Height 

11 

mm 

Width 

4 

mm 

Depth 

1 

mm 

Ia 

1395 

mm4 

ksw 

0.165e-5 

m 

Volt 

Nominal  CG  Location 

[5.5  0  0  ] 

mm 

Cmax 

4 

mm 

db 

3 

mm 

R 

15 

mm 

cqc^ 

V 

< 

[3.5  2  0  ] 

mm 

V 

< 

[  3.5  -2  0  ] 

mm 

Oiu 

45 

deg 

OtD 

45 

deg 

PTovpt  PrmiiFnrv  —  d1  —  /  weight 

113.61 

Hz 

mover  frequency  —  uQ  —  yj  5pCUa)lA 

Table  2.  Aerodynamic  forces  expressed  in  local  spar  and  body  frames. 

Figures  5  and  6  shows  the  commanded  and  actual  positions  and  attitudes  of  the  vehicle 
center-of-gravity  used  to  achieve  waypoint  following.  Tracking  performance  is  acceptable 
and  the  vehicle  is  capable  of  performing  the  desired  maneuvers.  Figures  7  shows  the  left  and 
right  wing  fundamental  and  split-cycle  frequencies  used  during  this  simulation.  There  are 
small  differences  between  the  fundamental  and  split-cycle  frequencies  of  the  left  and  right 
wing.  These  differences  result  in  a  roll  for  the  vehicle  and  are  used  to  align  the  vehicle’s 
heading  with  the  next  waypoint.  The  voltage  applied  to  the  bob-weight  actuator  and  the 
location  of  the  vehicle  center-of-gravity  in  the  z  body-axis  direction  are  shown  in  Figure  8. 
The  nominal  z  body-axis  center-of-gravity  location  is  0.  Hence,  it  can  be  seen  that  the 
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center-of-gravity  is  shifted  by  a  max  of  about  —0.05  mm.  This  corresponds  to  a  —0.21  mm 
shift  in  the  z-position  of  the  bob- weight.  The  bob- weight  was  constrained  to  move  no  more 
than  0.25  mm  so  this  is  nearly  the  maximum  z  direction  shift  of  the  center-of-gravity.  These 
results  show  that  the  bob-weight  is  active  and  is  used  to  generate  pitching  moments.  The 
trajectory  that  was  chosen  for  this  simulation  is  acceptable  and  the  vehicle  has  sufficient 
control  authority  to  perform  the  necessary  maneuvers. 


Figure  4.  Desired  trajectory. 


VIII.  Summary 

This  paper  has  described  the  split-cycle  formulation  and  has  presented  a  derivation  of 
all  split-cycle  related  parameters.  The  effects  of  a  bob-weight  upon  the  body-axis  forces 
and  moments  were  computed  and  then,  expressions  for  the  cycle-averaged  control  deriva¬ 
tives  associated  with  the  variations  in  the  fundamental  wing  beat  frequencies,  split-cycle 
parameters,  and  bob-weight  position  were  developed.  The  expressions  were  linearized  about 
a  hover  flight  condition  and  a  control  allocation  strategy  was  proposed.  Insight  into  the 
controllability  of  this  vehicle  was  provided  by  interpreting  the  structure  of  the  control  ef¬ 
fectiveness  matrix.  The  analysis  presented  shows  that  the  use  of  split-cycle  constant-period 
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Inertial  Positions  vs.  Time  (desired  and  actual) 
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Figure  5.  Inertial  positions. 

frequency  modulation,  with  independently  actuated  wings,  allows  one  to  manipulate  the 
x-body  and  z-body  axis  forces.  The  y-body  force  is  not  directly  controllable  using  the  split- 
cycle  approach  as  presented.  Direct  manipulation  of  the  y-body  force,  given  the  physical 
design  constraints  of  the  proposed  vehicle,  would  require  that  the  frequency  of  the  oscilla¬ 
tors  driving  the  wings  be  varied  more  than  once  per  cycle.  Rolling  and  yawing  moments 
can  also  be  independently  varied  using  the  split-cycle  technique.  A  pitching  moment  can 
also  be  generated;  however,  it  cannot  be  generated  independently  of  the  z-body  force  or  the 
rolling  moment.  Thus,  a  bob-weight  was  introduced  to  cancel  these  undesirable  changes 
in  the  pitching  moment  when  manipulating  the  split-cycle  parameters  to  generate  z-body 
forces  and  rolling  moments.  Given  the  ability  to  independently  manipulate  5  out  of  the  6 
cycle-averaged  body-axis  forces  and  moments,  untethered  controlled  flight  with  insect-like 
maneuverability  appears  to  be  feasible. 
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Figure  6.  Attitudes. 
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Delta  RW&LW  (rad/sec)  Omega  RW&LW  (rad/sec) 


Fundamental  Wing  Beat  Frequency  and  Split  Cycle  Parameter  vs.  Time 


Figure  7.  Fundamental  and  split-cycle  frequencies. 
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Bobweight  Actuator  Voltage,  Z  CG  Location,  and  Bobweight  Z  Position  vs.  Time 
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Figure  8.  Voltage  applied  to  bob-weight  actuator  and  center-of-gravity  shift  in  z  direction  due 
to  bob-weight  movement. 
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